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Abstract 

The dependence of the J/tp yield on the transverse energy Et in heavy ion 
collisions is considered within the statistical cc coalescence model. The model 
fits the NA50 data for Pb+Pb collisions at the CERN SPS even in the high-£ T 
region (Et ^ 100 GeV). Here -Ey-fluctuations and EWosses in the dimuon 
event sample naturally create the celebrated drop in the J/ip to Drell-Yan 
ratio. 
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Measurements of the J / if) suppression pattern in heavy ion collisions had been proposed 
as a diagnostic tool for the deconfmement of strongly interacting matter at the early stage 
of the reaction |1| (see also 0] for a modern review). However, it was found that the J /if) 
suppression in proton-nucleus and nucleus-nucleus collisions with light projectiles (up to 
S+U) can be explained by inelastic collisions with the nucleons of the incident nuclei (the 
so-called 'normal nuclear suppression') ||. In contrast, the NA50 experiments with a heavy 
projectile and a heavy target (Pb+Pb) revealed 'anomalous' J /if) suppression (§-§[, i.e. 
an essentially stronger drop of the J /if) over Drell-Yan ratio R than it could be expected 
from a simple extrapolation of the light projectile data. The dependence of R on the 
neutral transverse energy Et, which is used as an estimator of the collision centrality, shows 
a two step behavior. The onset || of the anomalous J /if) suppression (deviation of the 
data downwards from the 'normal nuclear suppression' curve) supposedly takes place at 
Et ~ 35 GeV. The second drop has been noted at Et ~ 100 GeV. The initial interpretation 
U attributed these two drops to the successive melting of the charmonium Xc states (which 
contribute about 45% to the measured J /if) number) and primary J/^'s in the quark-gluon 
plasma. 

Other interpretations are, however, possible. Because of relatively high statistical errors 
of the data at low E T , the behavior around E T ~ 40 GeV can be equally well interpreted as 
smooth decreasing of the measured ratio R(Et) due to charmonium suppression by comoving 
hadronsQ ||. Discussing the second drop of the J /if) to Drell-Yan ratio, one has to take 
into account that the neutral transverse energy Et estimates different quantities in the 
regions below and above Et ~ 100 GeV. At Et ^ 100 GeV, it is proportional with good 
accuracy to the number of nucleon participants N p and provides a reliable measure of the 
collision centrality. The large Et 'tail', however, is produced mainly due to fluctuations of 
the transverse energy at almost fixed (nearly maximal) centrality. The influence of these 
fluctuations on the J /if) suppression can partially explain the second drop. Better agreement 
of the co-mover model with the data in the region Et ^ 100 is achieved when a loss of 
transverse energy in the dimuon data sample with respect to the minimum bias one is also 



taken into account 10 



The Et fluctuations appear to be important in the deconfmement scenario as well. It 
was found Jll| that while the first drop can be explained by melting of x c 's, the onset of 
primary J / if) melting does not produce any visible structure on the curve. The second drop 
appears due to the fluctuations of Et- 

All mentioned approaches [|1],[7|-[TT| , despite of their differences, have one common feature: 
charmonia are assumed to be created exclusively at the initial stage of the reaction in primary 
nucleon collisions. Any subsequent interactions (with sweeping nucleons, co-moving hadrons 
or quark-gluon medium) may only destroy them. Despite of rather successful agreement with 
the J /if) data, such a scenario seems to be in trouble explaining the if)' yield. The recent 



lattice simulations |12| suggest that the temperature of if)' dissociation T d {if)') lies far below 
the deconfmement point T c : T<i(if)') ~ 0.1-0.2T C [13|]. Therefore, not only the quark-gluon 
plasma, but also a hadronic co-mover medium should completely eliminate if)' charmonia in 



1 Later it was found that even the sweeping nucleons of the colliding nuclei might, under certain 
assumptions, produce the observed suppression M. 
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central Pb+Pb collisions at SPS. However, the experiment revealed a sizable ip' yield (see, 
for instance, [14JL51 ). It was observed [16| that ip' to J/ip ratio decreases with centrality 



only in peripheral lead-lead collisions, but remains approximately constant and equal to its 
thermal value at sufficiently large number of participants N p > 100. 

A natural way to explain the ip' yield in nucleus-nucleus collisions without running into 
contradiction with the lattice data is to assume that a charmed quark and an antiquark that 
have been created at the initial ('hard') stage of the reaction can coalesce to form a hidden 
charm meson at a later stage, close to the chemical freeze-out point. This possibility has 
been considered in the literature |T7|-p2| within various models. 

The statistical cc coalescence model of J / ip production suggested in will be con- 



sidered in the present paper. The charmed quark- antiquark (cc) pairs are assumed to be 
created at the initial stage of the nucleus-nucleus reaction in hard parton collisions. Cre- 
ation of cc pairs after the hard initial stage as well as their possible annihilation is neglected. 
Therefore, the number of charmed quark-antiquark pairs remains approximately unchanged 
during the subsequent evolution. They are then distributed over open and hidden charm 
particles at the final stage of the reaction in accordance with laws of statistical mechanics. 
It is postulated that the deconfined medium prevents formation of charmonia at early stage 
of the reaction. In fact, it was found |23] that J/ip becomes unbound slightly above the 
deconfinement point T^J/i/j) « 1.1T C . Therefore, all charmonia observed in central Pb+Pb 
collisions supposed to be created at the final stage of the reaction. 



It was shown in our previous paper [24] that the statistical cc coalescence model is 



in excellent agreement with the NA50 data for (semi) central (N p > 100) Pb+Pb collisions. 
However, the analysis of pi] was restricted to the region Et ~ 100 GeV, where Et estimates 
the collision centrality. Fluctuations of Et were not taken into account. The aim of the 
present paper is to extend the applicability domain of the model to E T > 100 GeV and 
describe the second drop of the anomalous J/ip suppression. Two effects appear to be 
important in this region: 

• Influence of fluctuations of thermodynamic parameters of the system on the J/ip pro- 
duction. 

• A loss of neutral transverse energy, Et, in the event sample collected with the dimuon 
trigger relative to the minimum bias events. 

As far as charmed (anti)quarks are created exclusively at the initial 'hard' stage of the 
reaction in collisions of primary nucleons, the average number of cc pairs (cc)ab(&) created in 
a nucleus-nucleus (A+B) collision at fixed impact parameter b is proportional to the number 
of primary nucleon-nucleon collisions and can be found in Glauber's approach: 

(cc) AB(b) = ABag N T AB (b). (1) 

Here Tab(P) is the nuclear overlap function calculated assuming a Woods-Saxon distribution 
of nucleons in nuclei |25| (see e.g. Appendix of f24j). a^ N is the production cross section of 
a cc pair in a nucleon-nucleon collision. There are indirect experimental indications |26]| and 



a phenomenological model [27] suggesting an enhancement of the open charm in heavy ion 
collisions with respect to a direct extrapolation of nucleon-nucleon data. Therefore, a^ E N in 
nucleus-nucleus collisions is not necessarily equal to its value measured in a proton-proton 
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experiment. We consider free parameter. Its value is fixed by fitting the NA50 

data. 

The primary nucleon-nucleon collisions in a nucleus-nucleus reaction are independent 
and the probability to produce a cc pair in a single collision is small. Therefore, event- 
by-event fluctuations of the number of cc pairs in nucleus-nucleus collisions at fixed impact 
parameter b are Poisson distributed. Assuming exact conservation of this number one obtains 
the following statistical coalescence model result for the average number of produced J/ip's 



per A+B collision [19 



( J /WAB(b) = {CC)AB(b) [I + {cc)abQ>)\ (^(T /i B )/2) 



+ . . . 



(2) 



Where T and \i B are, respectively, the system temperature and baryonic chemical potential. 



The dots in (|2[) state for higher order terms with respect to the ratio 



which can 



be safely neglected at N p > 100. Here No is the total open charm multiplicity calculated 
within the grand canonical ensemble of the equilibrium hadron gas model: 



N o = y Yl n j (T,fj lB ) = Vn {T,iJ,B)- 

j=D,D,D*,D*,... 



(3) 



The sum in (||) runs over all known (anti)charmed particle species |£8]. The total J/ip 
multiplicity Nj°^(T, fi B ) includes the contribution of excited charmonium states decaying 
into J/ip: 



jvrtot 



v E Rti 

j=J/ip,xi,X2,i>' 



J/^)n {T^ B )^Vn%{T^ B ). 



(4) 



Here R(j — > J/if>) is the decay branching ratio of the charmonium j into J/if>: R(J/if) 
J I if)) = 1, R(xi -> J /ti)) » 0.27, R(x2 -> J I if>) ~ 0.14 and R{$' -> J/ijj) w 0.54. The 
volume parameter V is fixed by the condition of the baryon number conservation: 



N P (b) = V J2 bjTij(T, n B ) = Vn B (T, (a b ) 

j=N,N,AA,- 



(5) 



Here bj is the baryon number of particle type j and the sum in (||D runs over all known 



(anti)baryon species ||28|| . The number N p (b) of nucleon participants at fixed impact param- 
eter b is found in Glauber's approach (see e.g. Appendix of [p4|)p|. 
The particle densities in are found from: 



nj(T, n B ) 



dj 
2^2 



k dk 



exp 



/m| + k 2 — fij 
T 



± 1 



(6) 



where rrij and dj are, respectively, the mass and degeneracy factor of particle type j; fij = 
bjUB + SjUs+Cjfi'C is the chemical potential with Sj and Cj being, respectively, the strangeness 



2 In what follows we neglect fluctuations of the number of nucleon participants at fixed impact 
parameter. 
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and charm of particle type j. The strange and charm chemical potentials us and /ic are 
fixed by the requirement of zero net strangeness and charm in the system. 

As was already discussed above, the analysis JOJ of the recent lattice data |L2] has 
shown that ip' cannot exist in hot hadronic medium. To explain the observed ip' multiplicity 
one has to assume that charmonia are formed close to the chemical freeze-out point. (For 
(semi)central Pb+Pb collisions at SPS the chemical freeze-out (almost) coincides with the 
hadronization.) The thermodynamical parameters T and lib at chemical freeze-out can be 
estimated by fitting the hadron gas model to the hadron yield data in Pb+Pb collisions at 
SPS P^-pT|. This procedure, however, provides only average (in a sense that will be defined 
later) values T and ~p B . In fact, the thermodynamic parameters are not the same for all 
events, but rather they are subject to fluctuations of dynamical origin. 

Only a fraction of the initial energy of the incident nuclei is converted to thermal energy 
E t h (including the rest masses) of the produced hadron gas. The rest is redistributed to 
kinetic energy of the collective flow. The contribution of each nucleon participant to E t h is 
not fixed. It fluctuates in accord with a hitherto unknown dynamical mechanism. 

Suppose one is able to measure both N p and E th . Then, events can be grouped in bins 
with different thermal energy per participant, E t h/N p . Different thermal energies per net 
baryon number cannot be obtained at the same values of the temperature and baryonic 
chemical potential of the hadron gas. Therefore, fitting the hadron gas model to the hadron 
yield data in each bin, one would obtain different values of T and \i B . In reality, however, such 
a selection of events is not feasible, therefore no experimental information about fluctuations 
of T and \xb is available. The procedure of [^9H3T| has to mix events with different E t h/N p 
and estimates average values T and ~p B . One can, however, make a reasonable assumption 
about the fluctuations. It was found that at different initial energies of colliding nuclei, the 
produced hadron gas chemically freezes-out at constant total thermal energy per hadron 

H: 



E th (T, n B )/N h (T, h b ) = const « 1 GeV. (7) 

Different initial energies result in different thermal energy per participant. Therefore, it is 
natural to assume that the same criterion is valid when E t h/N p fluctuates at fixed initial 
energy. 

At given T and the thermal energy of the hadron gas can be found from 

E th {T if jL B ) = V J £e j (T, flj ), (8) 

j 

where Ej(T,fiB) is the contribution of the hadron species j to the thermal energy density: 

-l 



d ■ f°° i 

ej(T,fi B ) = ^2 / k 2 dk^m] + V 



, y m i + k2 ~ to , , , 
exp - — ± 1 



(9) 



The number of hadrons Nh(T, //b) is given by 

N h (T, f i B ) = V 1 £n j (T tflj ). (10) 
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The sum in equations (|8|) and fllOD runs over all known stable hadron and resonance species 
2|. 

Equation ([?[) defines the curve in the (T, fi B ) plane along which the freeze-out parameters 
fluctuate. To compare our model with the experiment the relation of each point on this curve 
to the transverse energy Et of neutral hadrons at given impact parameter b is needed. The 
point (T,~p B ) corresponds to the average value Et, which in the framework of wounded 
nucleon model is proportional to the number of participant nucleons: 

E T = qN p (b). (11) 

The fluctuations of Et are known to be Gaussian-distributed |[34|| : 

pnq 2 aN p (b) \ 2q 2 aN p {b) J 

The subscript 'MB' states for 'minimum bias'. This probability distribution is valid for the 
event sample collected with the 'minimum bias' trigger, which includes allQ events having 
the given value of Et- The parameter values q = 0.274 GeV and a = 1.27 [|35| are fixed 
from the minimum bias transverse energy distribution ||. The deviations of Et from its 
central value are mostly due to change of the number of neutral hadrons rather than due to 
change of average energy per hadron. Therefore, Et is approximately proportional to the 
total multiplicity of final hadrons: E T ~ iV£ ot (T, /i B ). Due to resonance decay, Nl ot (T, n B ) 
is larger than the number of particles in the hadron gas at freeze-out Nh(T, fi B ). Namely, 
every term in the sum ( |T0"D should be multiplied by a factor gj > 1: 

{1 for stable and weakly decaying particles; , > 

J2k kR{j —>■ kh) for resonances. 

Here R(j —>■ kh) is the probability that the resonance j produces k hadrons in the final 
state. 

( E ) 

Hence, the average number of J/ij) mesons, \J / AB{b)i produced in an A + B collision 
at fixed impact parameter b and transverse energy Et can be found in the following way: 

1. T and fi B are found as a solution of coupled transcendental equations 

E th (T,fx B )/N h (T,fi B ) = E th {T,JI B )/N h {T,J[ B ) , (14) 
Nl ot (T, /x B )/iV* ot (T, Ti B ) = Et/Et . (15) 

The first of the above equations expresses the freeze-out criterion (|7|): constant thermal 
energy per hadron. The second one appears due to the proportionality between the 
multiplicity of final particles and the measured neutral transverse energy. 

2. The values of T and \i B found from (|i~4D and (|15D are substituted into the formula (0) 
giving the desired value of (J /^)ab%y 



3 In contrast, the dimuon trigger requires production of a dimuon pair. 



6 



In real experiment, the impact parameter is not fixed. The relevant quantity, the differ- 
ential cross section of the J / ip production in collisions of nuclei A and B with respect to the 

( F ) 

transverse energy Et, is obtained by integrating {J/wIabIv) over ^ ne i m pact parameter b: 



dEn 



2vr r dbb (J/V)^) Pj/4E T \b). (16) 

J 



Here Pj/^(ET\b) is the probability to produce neutral transverse energy Et at fixed impact 
parameter b, provided that a J/tp particle is produced in the same event. 

Formula fllBf ) refers to the total number of produced J/ip% while the measured quantity 
is the number of pairs originating from their decay and satisfying certain kinematics 

conditions. Therefore, the cross-section ( fL6|) should be multiplied by the probability of J/ip 
decay into a dimuon pair B J J* = (5.88 ± 0.10)% O and by a factor 77, the fraction of the 
dimuons satisfying the kinematical conditions of the NA50 spectrometer. Because the con- 
sidered J/i/j production mechanism is completely different from the "standard" production 
in hard collisions, the rapidity distribution of J/^'s may be very different from that in N+N 
collisions. Consequently, r\ should not be necessarily equal to r\ N ^ « 0.24 which is obtained 
from Schuler's parameterization [[3(J of proton-proton data. Therefore, r\ will be treated as 
one more free parameter. 

Finally, to make a comparison with the NA50 data, the J/tp cross section should be 
divided by the Drell-Yan cross section. Similarly to ([16]), it is found from 



d(J AB- 



>DY' 



dEn 



J d 2 b (DY')AB( b ) P D Y(E T \b). (17) 



Here quantity Pdy (E T \b) is similar to Pj/^(E T \b), (DY') A B{b) is the average number of Drell- 
Yan pairs produced in an A+B collision at fixed impact parameter b. The prime means that 
only the pairs satisfying the kinematical conditions of the NA50 spectrometer are taken into 
account. Similarly to cc pairs, the number of Drell-Yan pairs is proportional to the number 
of primary nucleon-nucleon collisions: 

(DY')AB(b) = ABo™T AB (b), (18) 

where the production cross section of Drell-Yan pairs in nucleon-nucleon collisions. 

This quantity is isospin dependent, therefore the average value should be used: 

AB 

aDY ' = -AB' (19) 

For the case of Pb+Pb collisions, A = B = 208 and a^f, h = 1.49 ± 0.13 /ib [0. 
Hence, the following formula will be used to fit the NA50 data: 




(20) 



As the first step, we assume that both Pj/^(E T \b) and P DY (E T \b) in (|TB[) and ( |T7D 
coincide with the 'minimum bias' probability ([12]). The result of the calculations is shown 
in Fig. |I] with the dotted line. Calculations with T = T and /i# = ~p B , i.e. ignoring the 
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fluctuations, were done for comparison (the thin dashed line). As is seen, the influence of 
fluctuations becomes essential only at E? > 100. This can be explained as follows. If E? is 
not very large, Et « qN p (0), the integration in ([16]) runs over the region where Et/N p > q, 
as well as over the one with Et/N p < q. The central point of the Gaussian fll2|) , Et/N p = q, 
corresponds to T = T and Hb — ~Pb- Deviations of T and \ib in both directions from their 
central values almost cancel each other. Therefore, the result does not differ essentially from 
that obtained at T = T and \ib = V-b- 

In contrast, when Et approaches the value qN p (0) or exceeds it, the cancellation does 
not take place any more. At large Et, the region Et/N p < q falls into nonphysical domain 
N p > N p (0). Only the tail of the Gaussian (0) with E T /N p > q contributes to the integral 



(Pf). This corresponds to fluctuational tail with T > T and hb <~Pb- 

Although both the thermal 'density' of J /if), r^fj^iT, as well as that of the open 
charm, no(T, /j,b), are very sensitive to the temperature, their ratio in the formula @, 

Iff^T^ B ) 1 n « ( 



A%(T,M 1 n%(T,fM B ) 
(N (T,fi B )/2) 2 V{n {T^ B )/2Y' 



(21) 



is much less sensitive. The additional J/ip suppression takes place mostly due to increasing 
of the system volume V. Smaller \i& leads to smaller baryon density tib(T,/j,b) in @- At 
fixed N p this results in a larger volume. Decreasing the number of J /if} with increasing 
volume at fixed number of cc pairs is intuitively clear: the larger is the volume the smaller 
is the probability that c and c meet each other and form a quarkonium state. 

As is seen from figure [l], the suppression due to the T and \ib fluctuations is not strong 
enough to explain the experimental data at Et > 100 GeV. The reason is that we have 
assumed that Pj/^(E T \b) = PMB(E T \b). In fact, however, production of a J /if) particle costs 
a fraction of energy of the colliding nuclei. Therefore, the multiplicity of ordinary hadrons 
(and consequently Et) in a J /if} event is, in average, smaller than that in a 'minimum bias' 
event (i.e. without registering J/if>)f\ at the same impact parameter. Hence, as was pointed 
out in [JID| , the probability Pj/^(E T \b) is slightly different from P MB (E T \b). Similarly to 
| |10| |, we assume that the nucleon pair that have produced a J /if) meson does not contribute 
to the transverse energy. Therefore, the probability Pj/^(ET\b) is obtained from formula 
(|12|) by the replacement N p (b) — > N p (b) — 2. This induces a tiny (about 0.5%) shift of the 
Et probability distribution. 

The transverse energy distribution of the Drell-Yan event sample Pdy{Et\^) should also 
differ from Pmb(-B'tI^) because of the energy loss for the Drell-Yan pair production. However, 
in the large Et region, where the mentioned difference is important, our result is to be 
compared with the data obtained from the minimum bias analysis. This means that the so 
called "theoretical Drell-Yan" instead of the physical one was used as a reference for J/ if) 
suppression pattern. This "theoretical Drell-Yan" was obtained neglecting the Et loss in 



4 The 'minimum bias' probability distribution Pjv/b(-Et|^) is basically related to the events without 
J/ip or a Drell-Yan pair in the final state. Only a tiny fraction of Pb+Pb collisions results in 
production of charmonia or Drell-Yan pairs, therefore their influence on the shape of the minimum 
bias probability distribution is negligible. 
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the dimuon data sample. Therefore, to be consistent with the data presented by the NA50, 
we put P DY (E T \b) = P MB (E T \b). 

The result of the calculations taking into account the modification of Pj/ 7 j,(Er\b) (ignor- 
ing, for the moment, fluctuations of T and is shown in figure [l] with the dash-dotted line. 
Again, a sizable effect (compare with the thin dashed line) is observed only at Et > 100. 
This is explained by the following properties of the Gaussian distribution. Let us consider 
a ratio of two Gaussians. One of them is shifted with respect to another one: 

exp j- [x - (x - A)] 2 } 

1 J =exp{-2A[x-(x -A)]}. (22) 



cxp 



{- [x - x ] 2 } 



If the shift A is small, the ratio ( p2|) is very close to 1 near the central point xq. But, due 
to exponential behavior, it differs essentially from unity at large distances from xq, which 
corresponds to tails of the Gaussians. As was mentioned above, the shift of Pj/^Et^) with 
respect to PMB{Er\b) is about 0.5%. Therefore, when the main contribution to the integral 
(|16|) comes from the vicinity of the central value of the Gaussian (12), the effect is small. 



In contrast, at large Et, when the integration runs only over the tail of the Gaussian, the 
effect becomes essential. 

Taking into account both effects, T and /xg fluctuations and Et loss in the J/ip event 
sample, one finds excellent agreement with the fitted data including high Et {Et > 100 
GeV) region (the thick solid line in figure |IJ) . The free parameters a^ 5 N and rj were fixed by 
fitting the model to the NA50 dataQ. 

As was mentioned above, the statistical coalescence model is not expected to describe 
small systems. This can be seen from the if)' data [|18|| . In the framework of this model the 
multiplicity of if)' is given by formula (0) with the replacement Nj°L — > N^. Consequently, 
the if)' to J /if) ratio as a function of the centrality should be constant and equal to its 
thermal equilibrium value. The experimental data |L5] (see also a compilation in ||18|| ) are 



consistent with this picture only at rather large numbers of participants JT6[. Therefore our 
fit is restricted to the region 

E T > 27 GeV, (23) 

which corresponds to N p > 100. Note that all NA50 data except the two leftmost points 
from the 1996 standard analysis set lie in the region (|23|). 

To check the robustness of our model with respect to possible uncertainties we use 
two independent freeze-out parameter sets. The first one [^] has been obtained assum- 
ing strangeness and antistrangeness suppression by a factor 7 S : 

T = 158 MeV, jl B = 238 MeV 7, = 0.79. (24) 



In the second one |31[ the complete strangeness equilibration has been supposed: 



5 Only the complete (i.e. including T and \ib fluctuations and Et loss) model has been fitted to 
the data. The calculations shown by the thin dashed, dotted and dash-dotted lines in figure ghas 
been done with the parameter set fixed from the complete model fit. 
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T = 168 MeV, ]I B = 266 MeV 7, = 1. (25) 

The predicted effective cc production cross section as well as the fit quality appeared to be 
insensitive to the freeze-out parameters: 

a™ = 35.7 fib , X Vdof = 1.06 (26) 
for both sets ( p4|) and fl25|). The parameter 77 only slightly depends on T and ~p B : 

17 = 0.134 ±M (27) 

for the set (f24|), and 

77 = 0.126 ±8;811 (28) 

for the set ©. 

As it was already pointed out in our previous papers fH||2"4] |, the model predicts a rather 
large enhancement of the total charm (by a factor of about 3.5 within the rapidity window 
of the NA50 spectrometer). A direct measurement of the open charm in Pb+Pb collisions 
at the CERN SPS would allow for a test of this prediction. 

Extrapolation of the model to low Et region reveals strong disagreement with the ex- 
perimental data (about 3.5 standard deviations at the leftmost experimental point). This 
suggests that the charmonium production mechanism in central and peripheral Pb+Pb 
collisions may be of a very different nature. The 'standard' mechanism, formation of pre- 
resonant charmonium states in primary nucleon collisions and their subsequent suppression 
in nuclear medium, seems to dominate peripheral collisions. In contrast, the present sta- 
tistical cc coalescence model assumes the formation of a deconfined quark-gluon medium 
in central collisions. This medium destroys all primary cc bound states. The final state 
charmonia are formed at a later stage of the reaction from c and c quarks available in the 
medium. The distribution of c and c over charmonia and open charm hadrons follows from 
the laws of statistical mechanics. 

In conclusion, we have shown that the statistical cc coalescence model provides a quan- 
titative description of the 'anomalous J /if) suppression' in Pb+Pb collisions at CERN SPS. 
In particular, the 'second drop' of the J /if) suppression pattern, at Et ~ 100 GeV, can be 
naturally explained within the model as the result of two simultaneous effects: the fluc- 
tuations of the freeze-out thermodynamic parameters and transverse energy losses in the 
dimuon events with respect to the minimum bias ones. The model predicts a rather strong 
enhancement of the open charm. Therefore a direct measurement of the open charm would 
be a crucial test of the statistical cc coalescence model at the SPS. 



6 A different possibility has been considered in ||38fl . It has been assumed that both direct and 
statistical coalescence mechanisms of J /if) production are present at all energies and at all collision 
centralities. No open charm enhancement has been permitted in |3£|| . Therefore the contribution 
of statistical coalescence mechanism has been found to be small at the SPS. Then, the primordial 
J/if) } s that survived suppression in the quark-gluon plasma dominate even in the most central 
Pb+Pb collisions. Such an approach is, however, hardly able to describe the centrality dependence 
of the J /if) suppression pattern. 



10 



ACKNOWLEDGMENTS 



The authors are thankful to K.A. Bugaev, P. Bordalo, M. Gazdzicki and L. McLerran 
for comments and discussions. We acknowledge the financial support of the Alexander 
von Humboldt Foundation, DFG, GSI and BMBF, Germany. The research described in 
this publication was made possible in part by Award # UP1-2119 of the U.S. Civilian 
Research and Development Foundation for the Independent States of the Former Soviet 
Union (CRDF) and INTAS grant 00-00366 and INTAS grant 00-00366. 



11 



REFERENCES 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 



T. Matsui and H. Satz, Phys. Lett. B 178 (1986) 416. 
H. Satz, Rept. Prog. Phys. 63 (2000) 1511 ||hep-ph/000706^ 
M. C. Abreu et al, Phys. Lett. B 466 (1999) 408. 

M. C. Abreu et al. [NA50 Collaboration], Phys. Lett. B 410 (1997) 337. 
M. C. Abreu et al. [NA50 Collaboration], Phys. Lett. B 450 (1999) 456. 
M. C. Abreu et al. [NA50 Collaboration], Phys. Lett. B 477 (2000) 28. 
J. w. Qiu, J. P. Vary and X. f. Zhang, Nucl. Phys. A 698 (2002) 571 [|nucl-th/0106040 



C. Spieles, R. Vogt, L. Gerland, S. A. Bass, M. Bleicher, H. Stocker and W. Greiner, 
Phys. Rev. C 60 (1999) 054901 ||hep-ph/9902 , 337| ; 

J. Geiss, C. Greiner, E. L. Bratkovskaya, W. Cassing and U. Mosel, Phys. Lett. B 447 
(1999) 31 [ |nucl-th/9803008l ; 

N. Armesto, A. Capella and E. G. Ferreiro, Phys. Rev. C 59 (1999) 395 §epj 
ph/9807258|| ; 



D. E. Kahana and S. H. Kahana, Prog. Part. Nucl. Phys. 42 (1999) 269. 
A. Capella, E. G. Ferreiro and A. B. Kaidalov, Phys. Rev. Lett. 85 (2000) 2080 jhep- 
ph/000230q . 



[10] A. Capella, A. B. Kaidalov and D. Sousa, |nucl-th/010502~T| . 
[11] J. Blaizot, M. Dinh and J. Ollitrault, Phys. Rev. Lett. 85 (2000) 4012 [|nucl-th/0007020 



[12] F. Karsch, E. Laermann and A. Peikert, Nucl. Phys. B 605 (2001) 579 ||hep-lat/ 00 12023 



[13] S. Digal, P. Petreczky and H. Satz, Phys. Lett. B 514 (2001) 57 ||hep-ph/010523l |. 
[14] M. C. Abreu et al. [NA50 Collaboration], LIP-96-04 Invited talk at 26th International 

Symposium on Multiparticle Dynamics (ISMD 96), Faro, Portugal, 1-5 Sep 1996. 
[15] M. Gonin et al., PRINT-97-208 Presented at 3rd International Conference on Physics 

and Astrophysics of Quark Gluon Plasma (ICPAQGP 97), Jaipur, India, 17-21 Mar 

1997 



[16] H. Sorge, E. Shuryak and I. Zahed, Phys. Rev. Lett. 79 (1997) 2775 [[hep-ph/9705329 



[17] M. Gazdzicki and M. I. Gorenstein, Phys. Rev. Lett. 83 (1999) 4009 [|hep-ph/9905515 



[18] P. Braun-Munzinger and J. Stachel, Phys. Lett. B 490, (2000) 196 | frucl-th/0007059| . 
[19] M. I. Gorenstein, A. P. Kostyuk, H. Stocker and W. Greiner, Phys. Lett. B 509 (2001) 
277 |frep-ph/ 00 101481 ; 

J. Phys. G 27 (2001) L47 |hep-ph/ 0012015 . 
[20] P. Csizmadia and P. Levai, |hep-ph/0008T95| ; 

P. Levai, T. S. Biro, P. Csizmadia, T. Csorgo and J. Zimanyi, J. Phys. G G27 (2001) 
703 |frucl-th/ 001 1Q23| . 



[21] S. Kabana, J. Phys. G G27 (2001) 497 ||hep-ph/0010228| . 

[22] R. L. Thews, M. Schroedter and J. Rafelski, Phys. Rev. C 63 (2001) 054905 fKep- 
| ph/0007323| . 



[23] S. Digal, P. Petreczky and H. Satz, Phys. Rev. D 64 (2001) 094015 ||arXiv:hep- 
| ph/01060Tii 

[24] A. P. Kostyuk, M. I. Gorenstein, H. Stocker and W. Greiner, Phys. Lett. B 531 (2002) 
195 ||arXiv:hep-ph/0110269 



[25] C. W. De Jager, H. De Vries and C. De Vries, Atom. Data Nucl. Data Tabl. 14 (1974) 
479. 

M. C. Abreu et al. [NA38 and NA50 Collaborations], Eur. Phys. J. C 14 (2000) 443. 



12 



[27] A. P. Kostyuk, M. I. Gorenstein and W. Greiner, Phys. Lett. B 519 (2001) 207 jhcp- 



ph/0103057 



D. E. Groom et al. [Particle Data Group Collaboration], Eur. Phys. J. C 15 (2000) 1. 

F. Becattini, J. Cleymans, A. Keranen, E. Suhonen and K. Redlich, [hep-ph/0002267| . 

G. D. Yen and M. I. Gorenstein, Phys. Rev. C 59 (1999) 2788 ||nucl-th/98080T2f . 

P. Braun-Munzinger, I. Heppe and J. Stachel, Phys. Lett. B 465 (1999) 15 ||nucl- 
th/99030lU|| . 

J. Cleymans and K. Redlich, Phys. Rev. Lett. 81 (1998) 5284 ||nucl-th/9808030 |. 

A. Bialas, M. Bleszynski and W. Czyz, Nucl. Phys. B 111 (1976) 461. 

D. Kharzeev, C. Lourenco, M. Nardi and H. Satz, Z. Phys. C 74 (1997) 307 |Eepj 

ph/96122171] . 



B. Chaurand, quoted in |TT 



G. A. Schuler, CERN-TH-71 70-94, [hep-ph/9403387 . 

M. C. Abreu et al. [NA50 Collaboration], Phys. Lett. B 410 (1997) 327. 

L. Grandchamp and R. Rapp, Phys. Lett. B 523 (2001) 60 ||hep-ph/0 103124 



13 




20 40 60 80 100 120 

E T (GeV) 



FIG. 1. The dependence of the J/ip over Drell-Yan ratio R on the transverse energy. The points 
with error bars are the NA50 data. The lines correspond to different versions of the statistical cc 
coalescence model (see text for details) and to the normal nuclear suppression mechanism. The 
normal nuclear suppression curve is obtained at a a b s = 6.4 mb, where a a b s is the absorption cross 
section of preresonant charmonia by sweeping nucleons. The vertical line shows the applicability 
domain of the statistical cc coalescence model N p > 100. 
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